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Abstract 

Resolvent metrics are generalization of the resistance metric and 
jo ■ provide unified treatment of heat kernel estimates of sub-Gaussian 

£> , type under minimal conditions. 
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1 Introduction 

Heat propagation is not only interesting on its own, but reflects the very 
intrinsic structure of the space where it does take place. We gained such 
a new insight by Kigami's resistance metric [TD]. Unfortunately the use 
of resistance metric is applicable only on recurrent spaces. In this paper 
we eliminate this limitation and extend his notion to transient spaces, in 
particular for a class of transient graphs. 

In the course of study of heat propagation the analogy between results on 
continuous and discrete spaces is utilized (see e.g. [5]) and switching between 
them become a powerful tool of the studies. That is why we think that it is 
useful if we tackle, in the present paper, the technically less involved random 
walk case and return to the measure metric space version in a forthcoming 
paper. 

Kigami's work and several other papers inspire the following questions. 

For any given measure (Dirichlet) space is there a "good" metric in which: 

• the elliptic Harnack inequality holds 

• a parabolic Harnack inequality holds (in conjunction with heat kernel 
estimates)? 

The presented results contribute to the answer of these questions. 
We introduce the resolvent metric which is direct generalization of the 
resistance metric and we make the following observations. 

• Under the resolvent metric the volume doubling property implies the 
elliptic Harnack inequality. 

• Under the resolvent metric volume doubling property turns to be equiv- 
alent to the parabolic Harnack inequality and two-sided heat kernel 
estimates in a fully local sense (c.f. [E]). 

The main results are given in Theorem I5.lf5.2ll5.3l and 16.11 The paper 
concludes with examples. 



The main result of the paper can be summarized as follows. We consider 
a weighted graph (T,p) and a random walk on it. We assume that for all 
p xy > we have P (x, y) > po > uniformly. We construct the resolvent 
metric p and consider B p (x,r), balls in p, their volume V p (x,r) and define 
the scaling function F (x,r) = (r 2 V p (x,r)) for a well chosen m. Denote 
f(x, .) = F~ 1 (x, .) and p n (x,y) = p n (x,y) +p n+1 (x,y) the sum of the 
transition kernel. 

Definition 1.1 We define a set Wq of scaling functions, F : T x [0, oo] — y 
R + -.there is a C > such that for all x &T,r > 

F(x,2r) <c 
F(x,r) ~ 

Theorem 1.1 Volume doubling holds (V p G Wq) for p with respect to p if 
and only if there are C > c > 0,/3 > 1,5 > and an F e Wq such that for 
all x,y G T and n > 



Pn (x, y) < —, — 77 rr exp 

V p (x,f(x,n)) 

and for p (x, y) < 5f (x, n) 

Pn {x, y) > 



F(x,r)\ "- 1 



n 



V p (xJ(x,n)) 
hold. 

2 Basic definitions 

We consider (r, p) , weighted graph, T is a countable infinite set of vertexes 
and p xy = p yx > a symmetric weight. Edges are formed by the pairs 
for which p xy > 0. We assume that the graph is connected. These weights 
define a measure on vertices: 

V ( x ) = Yl Vx,y 

yer 

and on sets A G T 

p(A)=J2p(z). 



Due to the connectedness p (x) > for all x. It is natural to define the 
random walk on weighted graphs, which is a reversible Markov chain given 
by the one-step transition probabilities: 



P(x,y) 



\X (x) 



In what follows we always assume the condition (po): there is a constant p > 
such that for all x, y with p xy > 

P (x, y) > p 

holds. One can define the transition operator P on cq (T) functions Pf (x) = 
^2P (x,y) f (y) . The inner product for co(T,p) is defined by (/, g) = 

If p is a metric, balls are defined with respect to it by 

B p (x,r) = {y : p(x,y) < r} 

Denote B = B p (x, r) the connected component of B p (x, r) containing x. The 
volume of the connected part B is denote by V p (x, r) = p (B p (x, r)). 

Definition 2.1 We say that the volume doubling property, (VD) holds if 
there is a C p > constant such that for all x e T , r > 



<C n . 



V P (x, 2r) 

V p (x, r) 

3 The resolvent metric 

In case of recurrent spaces Kigami's observation is that the effective resistance 
between two vertices R (x, y) is metric. The existence of the resistance metric 
has a particular consequence that, for any / in the domain of the Dirichlet 
form £ 

\f(x)-f(y)\ 2 <R(x,y)£(f,f). (1) 

If the volume of balls Vr (x, r) with respect to the metric R satisfies the 
doubling condition the following important estimate holds 

R(x,B R (x,r)) x r. 

In particular 

R(x,B c R (x,r))>cr (2) 



while R (x, B R (x, r)) < r is evident. One may recognize that ([2]) holds only 
for recurrent weighted graphs. This is a nice particular situation which has 
been successfully utilized in several papers to obtain heat kernel estimates 
and stability results ( [6] , [10] [TT] ) . Almost the same proof which leads to ([2]) 
results the validity of the Einstein relation in the form: 

E x (T BR ( x ,r)) x rV R (x, r) 

and that the elliptic Harnack inequality follows from the bounded covering 
condition. Having all that the heat kernel estimates and the parabolic Har- 
nack inequality follows. 

The crucial observations fail in the transient case, first of all (|2j) obviously 
does not hold, since R(x,B R (x,r)) — > Ro > and the rate of convergence 
can be understand from the decay of R {Br (x, r) , B R (x, 2r)). 

In several previous works resolvents are used with success to analyze 
transient walks and diffusions. The simplest resolvent is the following: 

oo 
n=0 

It is clear that it is monotonically increasing in m and may be infinite if P n 
decays polynomially. In probabilistic terms one may consider this sum as 
the average visit time of y by the increasing family of independent walkers 
which has n m members at time n. In independent walkers we mean here 
that on a given site some new walkers "born" (according to the expansion 
of the family tree) and start independent walk. In what follows we need a 
modified version of the resolvent which provides nice correspondence to the 
power of the Laplace operator while it has basically the same propertyies. 
We fix an m G N which will be specified later and reserved as the parameter 
of the resolvent. 

In [9] we started the utilization of polyharmonic functions, Green function 
as well as Green operators (or resolvents). Now we follow this direction and 
find a new metric for non strongly recurrent graphs (weakly recurrent and 
transient) which posses nice features. 

Denote P the transition operator on l x (T). Pf (x) = Y^ y ~ x ^ ( x > v) f (v)- 

Definition 3.1 The Laplace operator is defined as A = P — I . The Dirichlet 
form corresponding to the Laplace operator is given by 



£i(f,g) = (- 


-A/,£) = 


= ((!- 


-P)f,g) = 


x,y 


-f{y)){g{x)- 


-g(y))n x ,y 



Definition 3.2 For A,BcT,Ar\B = $ we define the resistance 

R(A,B)=ini[S(fJ):f\ A = lJ\ B = 0]- 1 
/ec 

Let A C T, A A = P A - I, (A A ) m = {-l) m (I - P A ) m the ra-th iteration 
of the Laplace operator for m > 1 integer. If A = F we drop it from the 
notation. 

Let us recall, that A (A) = inf/^ ^ iM1 ' " " |2 - " f " 2 



~Jl7]p — • \\J ll - IM lU 2 (r» 
The domain of the Dirichlet form on T is defined by J- m = T A (£„ 
{/ G I2 (r, //) ,£ A (f,f) > 0}, where the bilinear form £^ is defined as 

rt(f,9) = ((l-P A ) m f,9) b w 

The quasi resolvent metric on A is defined as 

i^(x,y) = S up { l/( ^7 /(y)|2 : /(*) *f(y),fe T„ 

Note that i?^ is decreasing in A since £^ is increasing by definition, conse- 
quently R m 

R m (x,y) = S up| ^__ZM^ :/(*) ^ /(„),/ e T m \ . (3) 

is existing. 

That has the equivalent forms 

R m (x,y) = sup{\g(x) - g (y)\ 2 : <£ m (g,g) < 1} 
and 

J R m 1 (A, J B) = inf{^ m (/,/):/G J- m ,^ m (/,/)> 0,/U(x) = l,/| B = 0}. 
The former one can be seen using q = , ■* 

y/£mV,f) 

Lemma 3.1 1. For any f G J 7 (£ m ) 

\f(x)-f(y)\ 2 <R m (x,y)£ m (fJ). (4) 

2. If T is connected and A fl B = t/ien 

< i?" 1 (A, B) < 00. 



Proof. The statements follow from the definition. ■ 

Lemma 3.2 IfAcBcDcT then 

R m (A,B c )<R m (A,D c ) (5) 

Proof. By definition for D c C B c 

R^(AB C ) = mi{S m (fJ):f\ A (x) = lJ\ Bc = 0} 
> mf{£ m (f,f):f\ A (x) = l,f\ D c = 0} 

■ 

Lemma 3.3 R m (x,y) is a quasi metric: 
for any x,y G T, 

R m (x, y) = R m (y, x) , (6) 

R m (x,y) = if and only if x = y (7) 

R m (x,y) < 2(R m (x,z) + R m (z,y)). (8) 

Proof. The first statement ensured by the definition. For the second see the 
end of the proof of [12] Proposition 3.1. The weak triangular inequality can 
be see as follows: 

Rm(x,y) = sup {\g (x) - g(y)\ 2 :0 <£ m (g,g)<l} 

9 

< sup {2 \g (x) -g(z)\ 2 + 2 \g (z) - g (y)\ 2 : < S m (g, g) < 1 } 

9 

< sup{2\g(x)-g(z)\ 2 :0<S m (g,g)<l} 

9 

+ sup {2 \g (z) -g(y)\ 2 :0<S m (g, g) < 1 } 

9 

= 2(R m (x,z)+R m (z,y)) 

■ 

The next result of Mac'ias and Segovia is essential in our work. 

Theorem 3.1 (ITS^)If X is a non-empty set and d is a quasisymmetric with 
constant K : 

d (x, y) < K (d (x, z) + d (z, y)) 

then, there is a metric p, such that 

d p (x, y)-p (x, y) 
with p = nr — t? ■ 

1 l+log 2 k 



Corollary 3.1 There is a metric p and C > 1 > c > such that for all 
x,y e r 

cp 2 (x, y) < Rm (x, y) < Cp 2 (x, y) (9) 

Based on this theorem we define balls with respect to p: B p (x, r) = 
{y : p (x, y) < r}. As in the case of the resistance metric it can be that the 
balls are not connected. Let B = B(x,r) = B p (x,r) C B p (x,r) be the 
connected subset of B p (x, r) containing x. With the same slight abuse of 
notation we shall use Br for the sets (balls) with respect to the quasi-metric 
R m . One can immediately observe that (T,p,p) satisfies volume doubling if 
and only if (T,R m ,p) does. In addition if the bounded covering property 
holds with respect to one of R m or p it holds for the other as well and it 
follows from volume doubling (c.f.. [H]) . 

Definition 3.3 The graph F with metric a satisfies the bounded covering 
condition if there is an integer M > such that for all x G T , r > the ball 
B a (x, 2r) can be covered at most M balls of radius r. 

Definition 3.4 The m-resolvent is defined for an integer m > as follows. 
Let Q m (n) = f^"^ 1 ) , A cT finite set and for x,y G A 

G ?n ( x , y) = ^Qm iji) P n (x, y) 

the corresponding Green kernel is g^ (x,y) = —j-^-G^ (x,y). 

The Green operators G A defined as usual. It is worth to observe imme- 
diately, that for m = G^ = I A and for m = 1 G^ = G A , the usual Green 
operator. For infinite A the resolvent operator may be unbounded and the 
Green function is oo. For finite sets due to the transience of the Markov 
chain with Dirichlet boundary, these objects are well-defined. 

Lemma 3.4 The Dirichlet Green kernel g A (x, y) for finite A C V is a re- 
producing kernel with respect to S m . 

Proof. Let u G Jv^wUc = 



MaM = ((- A ) mG -^y 



u 



s -h u 



= u (x) 
The next corollary is immediate. 



Corollary 3.2 

£m (gt ( X > •) , #m 0> •)) = £m (a?, ») 

Lemma 3.5 The minimal value in the definition of R m (x,A c ) of £ m (f, /) 
zs ta£en &y # (y) = gA^ gt fa 2/) anrf 

iC(*,A c )=^(x,x) (10) 

Proof. Let h be an other function with h(x) = 1 H\a c = then d = h—g, h = 
d + g 

£ m (h, h) = £ m (g, g) + £ m (d, d) + 2£ m (d, g) 

but E m (d, d) > while £ m (d, g) = cd (x) = 0. ■ 

Lemma 3.6 Assume (I,/?) has the bounded covering property (or (VD) ), 
then 

„2 



R m (x,B c p (x,r)) > 



cr 



Proof. From §§§ and (jSD we have that there is a c > such for s = 

cr 2 ,B R (x,s) C B p (x,r) 

R m (x, B c p (x, r)) > i? m (x, ££ (x, s)) . 

If 

i? OT (a;, B C R {x, s)) > cs (11) 

we are ready. Now we prove ffTTl) following the steps of [12] • Let B = 
B R (x, s) ,y, z G -Br (x, s) and i? m (y, z) < As, A < 1. Let us fix a c and a 
z £ B with cis < R m (x, z) < s. We consider < q z (y) = 9 D^ x ' y l < 1 m- 
harmonic function on D = B\ {z} with q z (x) = l,q z (z) = 0. By definition 
and the reproducing property of the Green kernel 

1 1 

£ m (qz, q z ) = q z = 77777 r < 



R m (x,z) R m (x,z) 

where R m (x, z) denotes the resolvent metric within B, while 

\lz(y)\ 2 = \q z (y)-q z (z)\ 2 < R m (x,y)£ m (q,q) 
< CR m (y,z) < C(Xs) 2 1 
~ R m (x,z) ~ ( Cl s) 2 < 2 

if A = Ai is chosen enough small. Note that volume doubling implies bounded 
covering of B R (x, s) \B R (x, Cis) . Let B R (zi, Ais) the set of covering sets 
(via the covering with smaller B p balls: B p (zi,cr) C B R (zi.Xis) balls with 

9 



some extra increase of the covering number), i = 1,...,K, Denote q(y) = 
mi^ q Zi (y) and q = 2 (q - §) l^for)- We have that q(x) = 1 and g (y) = 
on 5^ (x, s). Finally we obtain that 



R£(x,B%(x,s))<£ m (q,q)<4^£ m (q z .q z )< . „ , r 

^-^ mini -Rm (z, ^jj 



4M 4M 
< . 

CiS 



Corollary 3.3 

R m (x,B c p (x,r))^cr 2 (12) 

Proof. The lower estimate was given above, the upper one is almost imme- 
diate. We chose S = Cr 2 so that B R (x, S) D B p (x, r) . Let y G dB R (x, S) 
and apply (jSJ) for {y} C £>£. (a;, S) = 



2 



R m (x, B R (x, S)) < R m (x, y) = S = Cr 



4 The tail distribution of the exit time 

This section contains two key results. One establishes an estimate similar to 
the Einstein relation, the other presents the estimate of the tail distribution 
of the exit time. The novelty in the approach is that in the lack of the 
usual Einstein relation all the arguments should be accommodated to the 
m-resolvent. 

For brevity we will use the following notations: 

E m (A\x) = E x (Q m+1 (T A ) \X = x) 

E m (A) = max^gA E m (A\x) 

E m (x,r) = E m (x,r) = K m (B p (x,r)\x). We will use the particular 
notation for m — 0, 

E p (x, r) = E (B p (x, r) \x) is the usual mean exit time. 



Lemma 4.1 For a set A C T, x G A, there is a Cq > 1 such that 
m ,™ x Em (A) Cn m 



CE m (A) E m (A) 



10 



Proof. 

T A < 2n + I(T A >n)T A oQ n , 
T% < 2 m ((2n) m + I(T A >n)T A n oQ n ), 

where Q n is the time shift operator. From the strong Markov property one 
obtains with C = 2^1 

E m {A) < C 2 n m + CE x (I(T A >n)E Xn (T^)) 

< C 2 n m + C¥ x (T A >n)E m (A). 

E m (A) Cn m m /rTi 

< = r — + F x (T A >n) 



CE m (A) ~ E m (A) 

and the statement follows. ■ 

Let us recall here that under (VD) the scaling function H (x, r) = 
r 2 V p (x, r) has nice regularity properties. 

Corollary 4.1 If(VD) holds then there is a c such that ifn = (^CQ 2 E m (x, r)) ™ 

P* (T Bp (x,r) >n)> c . (13) 

Here C$ is given by the Lemma \4-1\ 

Theorem 4.1 If (T, p) satisfies (VD) then, for B = B p (x,r) 

Px {TB p (x,r) < n) < Cexp(-ck m (x,n,r)) (14) 

where k = k m (x, n, r) > 1 is the maximal integer for which 

^<q min E m (B p (y,j)), (15) 

K y&B p (x,r) V P V k/ / 

where q is a small constant (to be specified later). 

Definition 4.1 Let us define (3 m as the smallest possible exponent for which 

*"'<*■*> < C (*Y\ (16) 



r 2 V p (x,r) \r J 

and observe that f|T6|) equivalent to (VD) . 

Remark 4.1 There are several further trancripts of (fl4"|) . In the simplest 
case if r 2 V p (x, r) x r 13 ,B = B p (x, r) one has 

P s (T B <n)<Cex P (-c^V" 1 ). (17) 



11 



Remark 4.2 From (|17|) one can see that the estimate is weaker as m in- 
creases. If a lower estimate of the same form and magnitude is aimed, m 
should be chosen as small as possible. However it should be recognized, that 
the increase of m not only increase the upper bound but the probability on the 
left hand side of (TT4")) . 

Proof of Theorem 14.11 The proof follows the old, nice idea of [I] (see 
also PP Lemma 3.14). The only modification is that we use the very rough 
estimate: 

k 

m 



1 B p (x,r) d- / Ji 



j=l 



where r% is the exit time of dB p (^, i) ,£,i = -^n_i an d k > 1 will be chosen 
later. From Lemma 14.11 we have that with t — ? 

P (r < t) < p + at m (18) 

where p 6 NL 1 — el and a = = 2 T - -, . Let ri be such that P(t <t) = 
(p + at" 1 ) A 1. The relation ( TT8"j) can be rewritten as 



P (r m < s) < p + as 

E (exp (-Ar m )) < E (exp (-Xr] m )) <p + a\-\ 
From that point the proof can be finished as in pQ . ■ 

4.1 The Einstein relation 

The relation between the mean exit time of a ball, its volume and resistance 
is regarded as a key tool to obtain heat kernel estimates. In this section we 
obtain the corresponding relation with respect to the distance p assuming 
only volume doubling and existence of the m > integer. More precisely 
we show the following statements. 



Theorem 4.2 7/(r,/i, p) satisfies (VD) then, satisfies (ER) : 

E p (x, 2r) x [R m (x, B c ) V p (x, 2r)] 1/m (19) 

with B = B p (x, 2r) , E p (x, r) = E m (B p (x, r) \x) . 

Theorem 14.21 will follow from the next statement and from the tail esti- 
mate (I26D of the exit time. 



12 



Theorem 4.3 If ft satisfies (VD) p then, (ER) m : 

E m (x, 2r) x R m (x, 2r) V p {x, 2r) (20) 

holds, where B = B p (x, r) . E m (x, r) = E m (B\x),R m (x, 2r) = R m (x, B c (x, 2r)) 

Let us recall, that K m (B\x) = E (Q m +i (T B ) \X = x) and Lemma 8.4 
from [hi] . 

The first lemma is elementary. 

Lemma 4.2 Lei B = B p (x,r),T = T B then 

E x (Q m +i (T)) x E x . (T m ) . 

Proof. Let T = T B ( x , r )- Assume that r is large enough to ensure T > 2m 
i.e. Bd (x, 2m + 1) C B p (x, r) and obtain 

> " T^ > <2m+l (T) > ^ r-L- > C > — . 



ml m\ m\ \ m J \2m J 

For small values the inequality follows by adjusting the constants. ■ 
Of course the statement holds for arbitrary finite set as well. 

Remark 4.3 From the definitions, the Theorem \4-l\ (VD) and (ER) it 
is immediate that 

s *jg&,M B >fad) (21) 



k + 1 y€B p (x,r) 

> cq min E m (B p (y,r)) k~^, (22) 

y£B p (x,r) 

(* + !)*.-' > cq E^Ml (23) 

k + 1 > c H^ ' (24) 



where B — B (x,r) , which yields 



P, (T B <n) < Cex P [-c( ^h ^ ) (25) 



P.,. (T Lt <n) < C'oxp(- ,-(^^ir" ' j (26) 



13 



Lemma 4.3 (Feynmann-Kacc formula, c.f.. J9$ or JXffij Let f be a function 
on T, A C T, A > satisfying 

Af - Xf = in B. 

Then for any x G A, u = (1 + A) ,T = Ta 

f (x) = E x [u T f (X T )] 

and for any m > 

( [J2 oo n+m Q m (n) P A ] /) (x) = E x (Q m+1 (T) W T +"7 (X T )) . (27) 

Corollary 4.2 // we choose f = 1,X = we have from fl27|) t/iat 



£ m (A|x) = E x (Q m+1 (T A )) = J2 G i (x, y) . (28) 

yeB 

Lemma 4.4 For any A C T, x G A 

E m (A\x)<g*(x,x)n(A). (29) 

Proof. The proof follows from (J2SJ). ■ 

Corollary 4.3 

^(Alx^i^x,^)/^) 

and in particular for x G T, r > 

E m {x,r) < Cr 2 V p (x,r) 



Both statements direct consequence of 
Proof, (of Theorem 14 .3p . The upper estimate is provided by Lemma [4.41 
for the lower estimate we apply the proof of Proposition 4.2 in [T2]. Denote 
B = B p (x, 2r). We start with flU): If / G J 7 (£ m ) 

|/(x)-/(y)| 2 <i? m (x ;1 /)£ m (/,/) 
in particular let g (z) = g^ (x, z) and z <E B p (x, Sr) then 

\g(x) -g(z) | 2 < R m (x,z)£ m (g,g). 
From reproducing property of ^ (x, z) we have that E m (g, g) = R m (x, 2r) = 

|# (x) - g {z)\ 2 < C5 2 r 2 = C5 2 < C5 2 g (x) 2 . 
14 



We can choose 5 such that C5 2 < 2, and we obtain from g(z) < g (x) , that 
for z € B p (x, Sr) 

9^(x,z) > -g*(x,x). 
Now we finish immediately using the definition and (VD) . 

E m {x,2r) = ^2g^{x,y)fi{y)> ^ g*(x,z)/i(z) 

yeB zeB p (x,Sr) 

> 2 9 ™ ( x > X ^ V p ( X ' 6 ^ 

> cR m (B p (x,r),B c p (x,2r))V p (x,2r), 

where the last step follows from (TTOj) and (VD) . m 



After the above preparations the proof of Theorem 14.21 is immediate from 
Theorem 14.31 and the next Lemma. 



Lemma 4.5 



E p (x, r) x E m (x, r 



l/m 



Proof. Let B = B p (x,r) ,T = Tb From the Jensen inequality we obtain 
that for m > 1 

E m (x, r) = E x (Q m+1 (T BH(s , r) )) x E x (T m ) > [E p (x, r)] m . 

For the opposite estimate denote E = E m (x, r) and 

2c (£) 1/m 

E p {x,r) = J^P(T>n)> J^ P ( T > n ) 

n „= Co ( J B) 1 /™ 

> c Q E l/m P (T B > 2c E 1/m ) . 
Now we use Theorem 14.11 in particular (T2T)j) 

P* (T < n) < C exp f -c (^^ ^ ) (30) 

P, (T > n) 
l-P x .(T<n) 

> l-Cexp|-c(^)^)>l/2 



15 



if we chose n m > H (x,r) and c such that logC — c ( -£- \ = 1/2 i.e. 

1 ( \Prn~ 1 

c ° = 2 ( lo c-i/2 ) ^ ne P ro °f * s complete. 



Proof of Theorem (14. 2 j) . The statement is immediate from Lemma 14.21 
14.51 and Theorem 14.31 



5 Heat kernel estimates 

In this section we show that (VD) implies the of off-diagonal upper and near 
diagonal lower estimates. The proofs are adaptation of the ones developed 
in the works [H],[I2] and [TU] all, on resistance forms in case of recurrent (or 
strongly recurrent) spaces, graphs. 

It is standard to deduce a diagonal lower estimate from (fT3l , (see also in 
PU Theorem 6.2 ). 

Proposition 5.1 If (VD) holds there is a c> such that for all n > 

Q 

p 2n (x, x) > — — — --, (DLE) 

V p (x,f(x,n)) 

where f (x, .) is the inverse of F (x, .) in the second variable. 



Remark 5.1 If we consider the example, V p (x,r) x r a then H(x,r) 
r 2 V p (x, r) x r a+2 , F (x, r) x r~^~ , / (x, n) x n^+2 



P2n (x, x) > en ( a+2 > . 

Thus 

^n m p n (x,x) = oo. 

Remark 5.2 One may object that mil— (a ? 2 ) ) > — ^ holds for all m, that 
seemingly contradicts to the initial argument, which indicates that m had to 
be chosen enough large. Let us notice, that if m is not enough large then the 
whole above argument is meaningless, the resolvent operator G m is bounded 
and we can not subtract the needed asymptotic information from it. Among 
others the notion R m (x, y) similarly to the usual resistance metric R (x, y) 
in the transient case is meaningless, furthermore the key observations (jl]) 
and (jl]) are not in our possession. 
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5.1 Estimates of higher order time derivatives of the 
heat kernel and relation to £ m 

A fairly simple but powerful method is developed in the mentioned works 
(see also [7]). The key observation is the following (see for the simplified 
proof [10] Theorem 10.4). Without any further assumption for any finite set 

Acr _ 

, , 2R(x,A) V2 

Vn {X, X) < h 



n jjl (A) 



where R(x, A) = sup yeA R(x, y). We have the following version of the state- 
ment. 

Proposition 5.2 There is a C > such that, for any finite set A <Z T and 
x G A, n > integer 

Vn [X, X) < C[ + 



n m fi (A) 

where R m (x,A) = sup y€A R rn (x,y). 

Before we prove the statement we show how one can obtain the diagonal 
upper bound from it. 

Theorem 5.1 Assume (VD) then 

C 

Vn (a;, x) < — — — --, (31) 

V p (xJ(x,n)) 

where f (x, n) is the inverse of F (x, r) in the second variable, furthermore 

C 
p n (x,y) < . (32) 

y/Vp (x, f (x, n)) V p (y, f (y, n)) 

Proof. Let A = B = B (x, r) and choose r to have — m if' - = tv^- — r and 

H \ > / n m Vp(x,r) 

n = \Cr 2 V p (x, r)~\ then 

C 

P2[Cr^V p (x,r)]^ m V X i X ) ^ V~(x~r) 

( w C 

V2n {X, X) < 



V p (x,f(x,n))' 

This shows the statement for even n, for odd n it can be seen together with 
(132]) as in [9]. ■ 
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Definition 5.1 Let us introduce the time differential operator and its itera- 
tions for k > 1 , n > 

{Df)n = fn+2 — fn 

(D k f) n = (D(D^f)) n 

Lemma 5.1 Let f n (y) = p n (x,y) then there is a C > such that for all 
k>0,n>0,x eT 

(-l) k f^(x)<-^-j:f 2n (x) (33) 

(2n) 



(fc) 

2« 



Proof. From the spectral decomposition of p2 n (x, x) we know that h 
(—1) f n (x) > for all fc > and the same implies that the map n — > 
(—1) f h 2n+2 — h 2n ] (x) non- decreasing. We show the statement by induc- 
tion using a slightly stronger statement. Assume it holds for all < i < k, 

{A[S t n\) 

where q = 2~ ( - l+2 ' ) , S k = S k _i + c k , and note that for i = the assumption 
holds. 





= 


",(fc-l) , (fe-l)l 
li An ' i 4n+2 J 






< - 

t 


1 4[c fe nJ 
l^j 2. [^-2, ^4n+2-2iJ 




< - 

< - 


1 |\(fc-i) ^(k- 1 )" 

4|c*nJ r4n-4 L c fe nj ft 2n+2 j 

1 fcC*" 1 ) , 





~ 4 [c fe nj ' 4 ™~ 4 L C ^J 
now by induction, if m = n — \c k n\ 

Ak) . J_L(fe-l) <- J; l X / \ 

h in < 4m ^ 4m < 4m (4 L5fc _ imJ) fc-i^- 4 L5*-i-J W 

Let us recall that f 2k (x) is non-increasing in k and find that 

Am- A [S k _ 1 m\ = A(n- \c k n\ ) - A [S k -i (n - \c k n\ )J 

> A{n- \c k n\) - A[S k -m\ 

> An-A([c k n\ + [S k ^in\) 

> An-A[(S k ^ + c k )n\ 
= An- A [S k n\ 
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which leads to the needed inequality. 

kin ~ 77"^ ^kUn-4[S k n\ (X) . 

(4hS fc nJ) 
Finally observing that Sk = J^j =0 2~^ +2 ) we have that Sk < \ and we obtain 

(BSD. ■ 



Lemma 5.2 Again, if f n (y) = p n (x, y), then 

£ m (fn, fn) = ("I)" 1 (D m f) 2n (x) = (-I)™ f£> (X) 



Proof. 

but 



£m{fnjn) = {{I ~ PT fnJn) 
((I-P) m fn,fn) = (EC-!)' (7) /«+*>/« 



(-I)" P m /) 2 n (*) 



Proof, (of Proposition 15.2)) Let A C T be a finite set and choose y* so that 

P2n 0, y*) ■ = minp 2n (x, y) 



zeA zeA 

< ^2 B (V)<1, 

zgr 
and 

follows. Let us denote f n (y) = p n (x,y). By elementary estimates we have 
that 



- 2"7Ti + Rm ( x > ^) ^ m (^ 2n ' f 2n ) 

\X {A) 

I Q 
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where in the last step Lemma 15.11 is used. Solving this for p 2n (x, x) we 
obtain 

— / — 2 \ V 2 

/ \ ^ n Rrn\X,A) / 2 R (x,A) \ 

»«(*,*) < C 1 — ^— + ^_ + c 2 — ^— j (34) 

*°(^%W- (35) 

The proof is finished by noting that p2n+i {x, x) < p 2n {x, x) . ■ 

5.2 The off-diagonal upper estimate 

The off-diagonal estimate can be easily obtained from the diagonal one. 
Theorem 5.2 Assume (p ) , (VD) and (DUE) F then 

C 

V{xJ{x,n)) 

C ( (F p (x,d(x,y))\^ 

< ——. — — exp — c ' ■ ' 



V p (x,f(x,n)) y V n 

The proof is word by word the same as for Theorem 8. 5 in fT^j or an alter- 
native proof is combination of Theorem 8.6 and 8.10 in [Tffl . 

5.3 Lower estimates 

The next task is to show the Near Diagonal Lower Estimate (NDLE) F : 
There are 5 and c > such that , for all x G T, r > 0, y G B (x, r) , n > if 
p(x,y) < 5f(x,n) then, 

Pn (x, y) > 



V p (x,f(x,n)) 
Theorem 5.3 If(T,fi) satisfies (VD) p and(DUE) F then, (NDLE) F holds. 

Proof. First we prove 

. . c 

P2n [X, y) > — — ; r-. 

V p {x,f{x,n)) 

for x,y G T satisfying d(x,y) = mod 2. Let us choose r such that 
n = F p (x, r) = [r 2 V p (x, r)] 1/m and denote f 2n (y) = P2n (x, y) , then 
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\hn (x) ~ hn (y)\ 2 < Rm (s, 2/) £m (/, /) 



By Lemma [5.21 we have that 

l/jJn («) - /2n (v)\ 2 < R m (x, y) S m (/ 2n , / 2 „) = # m (x, y) (-1)™ /]„ 



and by Lemma 15.11 and the diagonal upper estimate 



,2 1 



I hn 0) ~ hn (V) | < Rm (x, y) ~P2n (x, x) 



Rm (x, Sr 



[F p (x,r)\ 



mP2n \X) X) 



_, Sr 2 1 2 

- C ^T77 V^ 2 " x ' x - 7 p 2n (z, a;) 

r z Vp{x } r) 4 

if 5 is small enough. The above inequality means that 

P2n(x,y) > -p 2n {x,x) > 



2"™^~ J -V p (x,ry 
Finally 

P2n+1 (>,?/) = V (X, Z) H(z)p 2n (z,y) 

> P0P2n(z,y) 

c 
V p (x, r) ' 



6 Stability 

In this section we show that (VD) implies the parabolic Harnack inequality 
via the two-sided estimates. It is an interesting by-product that in our sce- 
narios the volume doubling property implies the elliptic Harnack inequality. 
It is already known that (UEp) and (NDLEp) are equivalent to the F- 
parabolic Harnack inequality (c.f.. [H]) where F is properly regular function 
and both imply (VD) . 

Definition 6.1 The function class W\ is defined as follows. F e W\ if there 
are (3 > ft > l,C > c> such that for all R > r > 0,x eT,y e B (x, R) , 

m%^2< c («y. (36) 

\r J F(y,r) \r J 
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Definition 6.2 We say that (PH) f) the parabolic Harnack inequality holds 
for a weighted graph (T, fi) with respect to a function F G W\ if there is a 
constant C > such that for any x G T, R, k > and any solution u > of 
the heat equation 

d n u = Am 

onV> = [k,k+F(x,R)]xB(x,2R), the following is true. On smaller cylinders 
defined by 

V- = [k + -F(x,R),k + -F(x,R)}xB(x,R), 
andV + = [k + jF{x,R),k + F(x,R))xB(x,R), 

and taking (n_,a;_) G T>~ , (n+,x+) G T> + , 

d(x-,x + ) < n + — n_, (37) 

i/ie inequality 

u n S x -) < C , w n+ (x + ) 

holds, where the short notation u n = u n + m„+i was used. 

Remark 6.1 At present it is not clear how the > 1 condition follows from 
the assumptions, while we expect it holds if fi is (VD) . 

Definition 6.3 We say that the elliptic Harnack inequality, (H) holds with 
respect to [i, and p if there is a C > such that for all x G T , r > if h is 
harmonic on B p (x, 2r) : 

(I-P)h 0) = forx eB 0, 2r) 

then 

max h(z) < C min h (y) . 

zdB(x,r) ydB(x,r) 

Theorem 6.1 Assume that (r,/i) satisfies (p ), then the following state- 
ments are equivalent. 

1. (VD) p , andF eW t 

2. (UE) F and (NDLE) F hold for an F G W u 

3. (PH) F holds for F G W\ with respect to p-balls. 

Proof. The proof can be reproduced from the one of [H] Theorem 12.1. ■ 
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Theorem 6.2 Assume that (T,/S) satisfies (p ), then (VD) and F G W\ 

implies (H) . 

Proof. It is evident that (PH) F implies (H) and from Theorem 16.11 we 
know that (VD) p and F eW x implies (PH) F . 

■ 

Definition 6.4 Two weighted graphs T with \i and V with p! are roughly 
isometric (or quasi isometric) with respect to the metrics d, d! (c.f. ' t 8, Defi- 
nition 5.9]) if there is a map (ft from T to V such that there are a, b, c, M > 
for which 

-d (x, y)-b<d'((ft (x) , (ft (y)) < ad (x, y) + b (38) 

for all x, y G V, 

d'^(T),y')<M (39) 

for all y' G T' and 

-H(x) <u'((/)(x))<cli(x) (40) 

c 

for all x G I\ 

Theorem 6.3 The F -parabolic Harnack inequality is rough isometry invari- 
ant with respect to p and p' . 

Proof. We know that (VD) p on (r, //) if and only if (VD) , on (r, //) . The 
equivalence of (VD) and (PH) F on both graphs are given by Theorem 16. II 
and the statement follows. ■ 

6.1 Comments 

Kigami in [10] constructed a metric which is quasisymmetric to RVr (x, d (x, y)) + 
RVr (x, d (y, x)) . This procedure can be applied to p 2 (x, y) V p (x, p (x, y)) + 
p 2 (x,y)V p (y, p(x,y)) as well. All the conditions are satisfied to obtain a 
new metric a which is quasi symmetric to p. We know that (VD) implies 
(VD) a and by Kigami's result [TU] we have 

i \ ^ c ( f<r( x ,y) 

Vn [x, y) < ,, , 1 y^ exp [ -c 

Pn (x, y) > 



V a (x,g~ 1 (n)) \ \ n 

c 



V a (x,g- 1 (n)Y 
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where g^ 1 (n) the inverse of g (r) = r a V a (x,r) and a is the exponent 
determined by construction of a from p. The lower estimate and parabolic 
Harnack inequality follows for a as well. It should be noted here, that with 
the introduction of the second new metric the dependence from x in the 
exponential term is eliminated and F (x, r) replaced by g (r). 

7 Examples 

There are several examples of fractals and fractal like graph possessing two- 
sided heat kernel estimates (see [2]) and examples on which direction depen- 
dence destroy it [8]. The major contribution of resistance metric in showing 
heat kernel estimates is emphasized in [5] and [TO]. The same arguments 
apply to our work. 

Example 7.1 The graphs in Example 4- [6], are strongly recurrent. There, 
/3 > a i.e. strong recurrence is assumed, that is not needed anymore. For 
instance the high dimensional graphical Sierpinski carpet, can be handled by 
our method. 

Kigami constructed Q a family of fractal structures in [TU] . The structures 
can be discretized and get the so called graphical Sierpinski gaskets, graphs. 
The typical structures are recurrent but it is easy to lift up them and get 
transient graphs. In |[2j Barlow has Proposition 5 as follows. 

Proposition 7.1 Let a > 1 and a graph (T, (x) e Q which satisfies V (x, r) ~ 
r a and E (x, r) ~ r 13 with respect to the shortest path metric, furthermore the 
graph is very strongly recurrent (see the definition there). Let A > 0. Then 
there is a weighted graph (r,/2j such that V (x, r) ~ r a + x }j U f, E(x,r) ~ r 13 
and satisfies the elliptic Harnack inequality. 

The graph (l\yu) is product of (r, /x) and an ultrametric space. The 

choice of A > can ensure that /3 < a + A. From Proposition 3 of [2] 
we know that such graphs are transient hence the resistance metric is not 
applicable, while the resolvent metric does if volume doubling holds for it. 
Polynomial volume growth in the resolvent metric follow from the asymptotic 
spherical symmetry of the Green kernel and from polynomial volume growth 
in the original metric. 
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